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The velocity distribution of spheres rolling on a slightly tilted rectangular two dimensional sur- 
face is obtained by high speed imaging. The particles are excited by periodic forcing of one of the 
side walls. Our data suggests that strongly non-Gaussian velocity distributions can occur in dilute 
granular materials even in the absence of significant density correlations or clustering. When the 
surface on which the particles roll is tilted further to introduce stronger gravitation, the collision fre- 
quency with the driving wall increases and the velocity component distributions approach Gaussian 
distributions of different widths. 
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Recent simulations and experiments have shown that 
clusters are formed due to inelastic collisions in excited 
granular matter jj]-[3). The formation of clusters was 
shown to lead to non-Gaussian velocity distributions. 
Here we report that even in the absence of cluster forma- 
tion or significant density correlations, dissipation leads 
to strongly non-Gaussian velocity distributions. 

The statistical description of dissipative systems such 
as granular flows are of fundamental interest. Because 
thermal energies are irrelevant, the concept of "granular 
temperature" has been introduced based on the veloc- 
ity fluctuations of particles. A dissipative kinetic theory 
has been developed to describe rapid granular flows 0. 
In this theory, velocity distributions are assumed to be 
locally Gaussian. Deviations from Maxwell-Boltzmann 
are explained by averaging over local Gaussian velocity 
distributions with different widths which are inversely re- 
lated to the local clustering || . Given the success of this 
local Gaussian distribution assumption, we might con- 
clude that the velocity distributions can be considered to 
be approximately Gaussian for most practical situations 
where clustering does not occur. 

To test this idea, we study the velocity distributions 
of steel spheres rolling inside a two-dimensional box for 
which energy is supplied to the particles by oscillating 
one of the side walls. A similar system has been used 
previously to demonstrate the formation of clusters at 
high densities due to inelastic collisions Q. However, 
at low densities corresponding to longer mean free paths 
and lower collision rates, the particles appear to be ran- 
domly distributed and no clustering is observed. We use 
high speed and high resolution imaging to obtain accu- 
rate velocity distributions. Because one of the sides of 
the rectangular geometry is used to input energy and the 
particle-particle and particle-wall interactions are dissi- 
pative, the velocity distribution is asymmetrical. An ob- 
vious consequence is that the velocity fluctuations in the 
direction parallel to the driving wall motion is greater 
than in the perpendicular direction. We find that the ve- 
locity distribution is strongly non-Gaussian even in the 



direction perpendicular to the motion of the oscillating 
wall. By varying the frequency / and the angle of inclina- 
tion of the surface, we show that non-Gaussian velocity 
distributions can occur when the collision rate due to dis- 
sipative interactions is greater than the particle-driving 
wall collision rate. 

The experimental setup consists of one hundred 
0.32 cm diameter steel (or brass) spheres rolling on an 
optically smooth glass surface (29.5cm x24.7cm) which 
is enclosed by steel side walls. The setup is similar to 
that in Ref. |Q with a few enhancements. The coefficient 
of rolling friction of the spheres on the glass plate is less 
than 2 x 10 -3 . In contrast, the coefficient of sliding fric- 
tion is ~ 0.2. The particles are visualized using a SR-1000 
Kodak digital camera capable of taking 1000 frames per 
second (fps). The camera and lights are placed such that 
the centers of the spheres appear bright. The best pixel 
resolution of 512 x 480 can be obtained at 250 fps, which 
was used for most of the data presented here. A typical 
image is shown in Fig. [l]a. The particles correspond to 
~ 4 pixels in diameter and a centroid technique was used 
to find the position of the particle to within 0.01 cm. The 
particles with the highest velocities were determined by 
using images separated by 4 ms, and the velocity of the 
slowest particles were determined by using images sepa- 
rated by at least 12 ms. The positions of the slow moving 
particles do not change appreciably in comparison to the 
pixel size if a very high frame rate is used, resulting in 
substantial roundoff errors. We note that a combination 
of both high resolution and high frame rate is required 
to obtain accurate velocities. 

The driving wall is connected to a linear solenoid which 
is controlled by a computer, and the resulting motion of 
the piston is measured with a displacement sensor. The 
displacement of the piston as a function of time, d(t), for 
three typical driving frequencies is plotted in Fig. [l]b. A 
square waveform of fixed amplitude (4 volt) and various 
frequencies is used to obtain the displacements. 

By averaging over at least 10 5 images, we obtain veloc- 
ity distributions of the particles rolling inside the rectan- 
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gular box. The distributions of the velocity components, 
P{v x ) and P(v y ), are plotted in Fig. ||a and ||b respec- 
tively. The directions x and y and the origin are indi- 
cated in Fig. 0a. The surface is slightly tilted {9 = 0.1°) 
to continously excite the particles. Particles that are less 
than 2 cm from the side walls are excluded to avoid direct 
effects of the boundaries. We checked that if a narrow 
region in y is used, the form of the velocity distributions 
is independent of y. The plots in Fig. correspond to 
various driving frequencies. Note that the form of P(v x ) 
and P{v y ) does not depend strongly on the driving fre- 
quency /. The distribution P(v x ) is symmetrical about 
v x = 0, but is strongly non-Gaussian as can be seen from 
the poor fit to a Gaussian which is also shown in Fig. ||a. 

The distribution P(v y ) shown in Fig. |^b is non- 
Gaussian and asymmetrical. The reason for the asymme- 
try is not difficult to understand. The energy is supplied 
to the system by the driving wall. Therefore, particles 
move with higher velocities in the negative-?/ direction 
(indicated by the bump at v y ~ — 35cm/s). The second 
more subtle reason for the asymmetry is that the particles 
suffer an inelastic collision either with the wall opposite 
to the driving wall or with other particles before mov- 
ing in the positive y direction. Therefore, particles on 
average move with lower velocities in the +y direction 
than in the —y direction. Consequently, the particles 
also spend a longer time interval traveling in the +y di- 
rection than in the —y direction. Thus, asymmetrical v y 
distributions are obtained by averaging over time. Such 
effects were also observed in the distributions obtained 
numerically by Luding et al. ||. The form of P(v x ) is 
relatively independent of the energy input, and therefore 
its non-Gaussian behavior is surprising. 

In Fig. ||, we plot P(n), the distributions of the num- 
ber of particles n, in a 3 cm x 3 cm area averaged over 
the box and over all images. For comparison, the dis- 
tribution of randomly distributed particles is also shown 
(dashed curve). Our data for P(n) is similar to the ran- 
domly distributed case with a small deviation at higher 
n. Because of the dilute nature of our experiments, no 
effects of clustering are observed as reported in Ref. Q . 

The effect of gravity and the collision rate with the 
driving wall is studied by lifting the surface on which the 
particles roll by the tilt angle 9. The v x and v y distribu- 
tions are plotted in Fig. corresponding to four values 
of 9. The density distribution P(n) is no longer uni- 
form in y as in the 9 — 0.1° case, but decreases strongly 
with distance from the driving wall. (At the highest an- 
gle 9 = 6.7°, the particles do not reach the top wall.) 
The velocity distributions correspond to particles in a 
2 cm wide region at a distance y = 4 cm from the driv- 
ing wall to avoid the effects of the increase in potential 
energy and the decrease in the average v y at higher 9. 
The distributions P(v x ) and P(v y ) grow broader as 9 
increases. A Gaussian fit to the data corresponding to 
9 = 6.7° is shown in Fig. ||. The data for P(v x ) is rea- 
sonably described by a Gaussian especially at low v x . 
However, P(v y ) deviates systematically from Gaussian 



because of the asymmetry. Furthermore, the width of 
P(v y ) is large compared to the width of P(v x ) because 
energy is mostly supplied to the x component in an inelas- 
tic collision consistent with earlier observations of verti- 
cally vibrated granular matter |n| . We further note that 
P(v x ) averaged over all y is non-Gaussian for 9 = 6.7° 
because of the slight increase in the widths of the Gaus- 
sians as a function of y. This might also be the reason 
for deviations from Gaussian observed in velocity distri- 
butions of vertically vibrated monolayer layer of particles 

As 9 increases, the acceleration in the y direction in- 
creases. Therefore the particles fall back to the driv- 
ing wall faster as 9 is increased. This results in higher 
collision rates of particles with the driving wall leading 
to higher average velocities. Thus as the collision rate 
with the driving wall increases, the distributions become 
broader and more Gaussian. Furthermore, the total num- 
ber of inelastic collisions the particles undergo before re- 
turning to the driving wall decreases because fewer par- 
ticles hit the top wall at higher 9. 

We note that there is no discrepancy with the observa- 
tion of the broadening of the velocity distributions with 
an increase in 9 and the data plotted in Fig. ^. It might 
appear that changing the frequency of the driving piston 
should also increase the collision rate with the bottom 
wall and hence have the same effect. However, the distri- 
butions plotted in Fig. || do not depend strongly on the 
driving frequency. The lack of the strong dependence of 
P(v x ) and P(v y ) on / in Fig. || can be explained by the 
fact that fewer particles are struck by the forward moving 
wall at higher /. Thus the net energy input (and colli- 
sion rate with the bottom wall) appears to not depend 
strongly on /. 

Next we discuss the effect of rolling on the distribu- 
tions. In a previous study, it was noted that rolling ap- 
pears to lead to a lower effective coefficient of restitution 
r during collisions Q. Some theoretical work has been 
done studying rolling and collisions In Fig. 5a, we 
plot the r as a function of scattering angle of a rolling 
steel sphere on a glass surface colliding with steel bound- 
ary walls. The coefficient of restitution is calculated as 
the square root of the ratio of the final to the initial ki- 
netic energy. The data is distributed over a range of r 
which becomes smaller for higher <p. The large fluctua- 
tions are not related to measurement noise nor can they 
be accounted for by variations in the initial velocity (see 
Fig. 5b.) 

The mean value of r depends on <f> and is lower for 
higher <fi (indicated by the dashed line in Fig. 5a.) The 
mean value of r during a wall collision averaged over all <f> 
was 0.5±0.02. In contrast the average loss of energy dur- 
ing a particle-particle collision was lower and the mean 
r was found to be 0.65 by averaging over 20 collisions. 
In comparison, the coefficient of restitution of the steel 
sphere bouncing off a steel block after being dropped is 
0.93. The larger losses in kinetic energy occur mostly be- 
cause the particles slide as their angular velocity changes 
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to match their new translational velocity. The change in 
angular velocity was observed to occur within a centime- 
ter even at the highest velocities observed. It may be 
possible to approximately model the complications due 
to rolling by simply assuming a lower value of r. How- 
ever, an event-driven simulation using only a loss of ki- 
netic energy in the normal direction does not produce 
non-Gaussian distributions in v x , but yields asymmetric 
distributions for v y The subtle effects introduced 

by rolling may help us understand the effect of higher 
dissipation on velocity distributions. 

Limited experiments were also performed with brass 
spheres. The distributions obtained were identical to 
those corresponding to steel spheres under the same con- 
ditions. Thus the velocity distributions obtained are not 
very sensitive to the details of the surface properties of 
the spheres. 

The distribution of the speed of the particles is a com- 
bination of the non-Gaussian or approximately Gaus- 
sian velocity distributions of different widths (such as in 
Figs. ||andQ) and deviates from the Maxwell-Boltzmann 
form. Subtle effects on the velocity distributions such as 
the scaling of the temperature (the second moment of 
the distributions) in the near elastic limit has been stud- 
ied theoretically by Kumaran 0|, Grossman et al. [Q, 
and experimentally by Warr et al. jio]]. However, the 
deviations observed in our experiments are stronger as 
the loss of energy during a collision is significant com- 
pared to the energy of the particle. Therefore scal- 
ing around the Maxwell-Boltzmann distributions as dis- 
cussed in Refs. [ pii| , p^[p^ ] for r < I cannot be applied to 
our data. 

In conclusion, our experiments on excited dissipative 
granular matter rolling and bouncing inside a rectangular 
box have yielded velocity distributions that are strongly 
non-Gaussian in the absence of significant clustering and 
density correlations. Our data suggests that caution has 
to be used in assuming Gaussian velocity distribution and 
Maxwell-Boltzmann distributions even for rapid granular 
flows where position correlations are negligible. 
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FIG. 1. (a) Image of 100 steel spheres rolling on a glass 
surface. The surface is tilted at angle 9 to the horizontal 
(6* = 0.1°). (b) Example of driving wall displacement d(t) 
measured using a position sensor for frequencies / = 2 Hz, 
10 Hz, and 16 Hz. The displacements for the three cases are 
offset by 1 cm for clarity. The wall is driven using a linear 
solenoid. 



FIG. 2. The probability distribution of velocity compo- 
nents v x and v y [6 = 0.1°). The velocity distributions are 
strongly non-Gaussian. Note the peak at v y — —35 cms -2 
due to the driving wall. The corresponding driving wall dis- 
placement d(t) versus time is shown in Fig. nib. 



FIG. 3. The distribution of particles P(n) in a 3 cm x 3 cm 
area. The distribution corresponding to randomly distributed 
particles is also plotted. No significant correlation in the dis- 
tribution is observed. 



FIG. 4. (a) and (b) Velocity component distributions as a 
function of inclination angle 9 (/ = 10 Hz). All particles are 
in a 2 cm wide region at a distance y — 16 cm from the top 
wall. The v x and v y distributions approach Gaussians as 9 is 
increased. The curves are a Gaussian fit for 9 — 6.7°. 
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FIG. 5. (a) Effective coefficient of restitution r versus scat- 
tering angle <j> for a rolling particle colliding off the side walls. 
The dashed line is a least squares fit to the data, (b) The 
plot of r versus initial velocity v for 170 < (j> < 180 fluctuates 
significantly. 
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FIG. 1. (a) Image of 100 steel spheres rolling on a glass surface. The surface is tilted at angle 
9 to the horizontal (9 = 0.1°). (b) Example of driving wall displacement d(t) measured using a 
position sensor for frequencies / = 2 Hz, 10 Hz, and 16 Hz. The displacements for the three cases 
are offset by 1cm for clarity. The wall is driven using a linear solenoid. 
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FIG. 2. The probability distribution of velocity components v x and (9 = 0.1°). The velocity 
distributions are strongly non-Gaussian. Note the peak at v y = —35 cms -2 due to the driving wall. 
The corresponding driving wall displacement d(t) versus time is shown in Fig. lb. 
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FIG. 3. The distribution of particles P(n) in a 3 cm X 3 cm area. The distribution correspond- 
g to randomly distributed particles is also plotted. No significant correlation in the distribution 
observed. 
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FIG. 4. (a) and (b) Velocity component distributions as a function of inclination angle 9 ( f = 10 
Hz). All particles are in a 2cm wide region at a distance y = 16cm from the top wall. The v x and 
v y distributions approach Gaussians as 9 is increased. The curves are a Gaussian fit for 9 = 6.7°. 
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FIG. 5. (a) Effective coefficient of restitution r versus scattering angle (f> for a rolling particle 
colliding off the side walls. The dashed line is a least squares fit to the data, (b) The plot of r 
versus initial velocity v for 170 < (f> < 180 fluctuates significantly. 
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